Quasiparticles of spatially anisotropic triangular antiferromagnets in a magnetic field 
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The spectral properties of the spin- 1/2 Heisenberg antiferromagnet on an anisotropic triangular 
lattice in a magnetic field are investigated using a weak-interchain-coupling approach combined with 
exact solutions of a chain. Dominant modes induced by interchain interactions in a magnetic field 
behave as quasiparticles which show distinctive features such as anomalous incommensurate ordering 
and high-energy modes. In terms of them, various unusual features observed in the anisotropic 
triangular antiferromagnet CS2CUCI4 in a magnetic field are quantitatively explained in a unified 
manner. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 



In the context of two-dimensional (2D) spin liquids 
[H, identification of gapless quasiparticles (QPs) emerg- 
ing in the absence of classical long-range orders in 2D 
frustrated magnets has been one of the central subjects 
in modern condensed-matter physics. Various possibil- 
ities of QPs different from conventional magnons, in- 
cluding those with fractional quantum numbers called 
spmons , have been discussed for 2D spin liquids 

in the absence or presence [3] of a magnetic field. Related 
to this issue, the anisotropic triangular antiferromagnet 
CS2CUCI4 has attracted a lot of attention especially be- 
cause of its spinon-like behavior observed in a zero field 
Meanwhile, it shows a host of intriguing features in 
a magnetic field. Among them are novel incommensurate 
(IC) ordering and excitation spectra which are hard 
to understand in terms of magnons in linear spin-wave 
theory Nevertheless, the conventional 2D magnon 
appears above the saturation field Hs 0. Obviously, 
QPs in the one-dimensional (ID) Heisenberg chain in a 
magnetic field 0, H, @] cannot consistently explain the 
above features since they show neither IC ordering nor 
2D dispersion relations by themselves. 

To explain the puzzling behavior in a magnetic field, 
we consider QPs which are defined as collective modes 
induced by interchain exchange processes and can be re- 
garded as bound states (BS) or anti-bound states (ABS) 
of ID QPs in a magnetic field. In zero field, related col- 
lective modes [lol, IHI, [13, may appear more or less 
similar to conventional magnons because of their nearly 
sinusoidal dispersion relations with gapless points close 
to those of magnons. However, the QPs in a magnetic 
field we consider in this Letter show distinctive features 
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FIG. 1: (a) Anisotropic triangular lattice, (b) Momentum 
regions of J'(fc)>0 (light blue) and J'(fc)<0 (light yellow) for 
J'>0. Solid circles denote ordering momenta for m=l/16, 
which shift towards the center of each light yellow plaquette 
as the magnetic field increases as indicated by arrows, (c)-(e) 
Comparisons with experimental results on CS2CUCI4. Solid 
lines are present results, (c) Magnetization curve. The dotted 
line denotes experimental results for _HJ|c in Ref. [l^ with g- 
factor 5=2.30 ^ and J=0.374 meV [ef. S^ax denotes 5^ at 
Hs ■ (d) Dispersion relation above Hs along dotted lines in (b) . 
Symbols are experimental results in Ref. 6;]. (e) Momenta of 
IC ordering. Symbols are experimental results at fcy = for 
_ff||cinRef. |J|, replotted using the magnetization curve in (c) . 



in a magnetic field 0, [^, 0, [3| are consistently explained 
as properties of such QPs. 

In this Letter, we consider the spin-1/2 antiferromag- 
netic Heisenberg model on an anisotropic triangular lat- 
tice in magnetic field H . The Hamiltonian is defined as 



which differ from those of conventional magnons [14| and 
those of ID QPs 0, i, H not only conceptually but also 
in appearance. In particular, we mainly focus our at- 
tention on the following features of the QPs: (1) mul- 
tiparticle crossover in a magnetic field, (2) frustration- 
induced IC ordering whose momentum strongly depends 
on the magnetization, and (3) excitation spectra with 
anomalous high-energy modes. We find, rather remark- 
ably, that various unusual features observed in CS2CUCI4 
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where Sr is the spin-1/2 operator at site r, and 5^ = 
Sr ■ The intra- and interchain couplings and unit 
vectors are denoted by J, J', x, and y as in Fig. [Ha). 
We assume J>\J'\. In particular, we take J'/J=0.34 
corresponding to that of CS2CUCI4 (J=0.374 meV and 
J'=0.128 meV W|). We denote the magnetization per 
site by m. The dynamical structure factors (DSFs) 
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are defined as S°'°'{k,uj)=^^M""{k,ei)S{uj—ei), where 
Af""(fe,ej) = |(fe,ej|S'^|GS)p for a=+, -, and z. Here, 
GS) and |fc, e,j) eacfi denote tlie ground state and an ex- 
cited state with excitation energy and momentum k in 
the anisotropic-2D system in a magnetic field. We also 
define S{k, u)=[S-+{k, uj)+S+-{k, uj)+AS''{k, w)]/6. 

We apply two weak-interchain-coupling techniques, us- 
ing exact solutions of the Heisenberg chain [16|. One is a 



random-phase approximation (RPA) [ll|, [ij, [13| : DSFs 
are obtained from dynamical susceptibilities approxi- 

mated as x''''ik,u;)^xfSikx,uj)/[l+J'ik)xfSik^,u;)/U 
with ^z = l and ^±=2, where the Fourier component 

of interchain couplings is J'(A;)=4J'cos -^cos -j- for 
anisotropic triangular lattices, and XiSi^xy^) are calcu- 
lated from DSFs in ID as lmxfS{K,Lu)=TTS^(k^^ for 
ti;>0. The S"^{kx,uj) can be calculated |g.lr7[[la .[ig| us- 
ing rapidities of Bethe-ansatz solutions [16| . It is known 
that spectral properties of the Heisenberg chain in a mag- 
netic field are mainly characterized by QPs called psinon 
(■0) and antipsinon (ip*) M, [a. m and a QP for a 2-string 
(a) in the Bethe ansatz [9|, We use excitations of 

up to 2tp24!* 0, 01, 0{L^) states of 2-string solutions 
and those of 3-string solutions Q in a chain with length 
L=320 for the RPA. These excitations occupy over 80% 
of the total spectral weight in L=320 except at very low 
fields [a]. 

The other technique is the one developed in Ref. [lo| : 
in the restricted Hilbert space spanned by exact eigen- 
states of chains, an effective Hamiltonian is derived as 

(1) 

where A°'{kx, ei) = y/ AI"^{kx, e^). It is known that DSFs 
calculated using the eigenstates of Eq. |T]) are almost the 
same as those of the RPA ^0|] . We discuss physical pic- 
tures mainly based on this formulation. We use dynami- 
cally dominant excitations of O(L^) states in an L=2240 
chain 0, Q and assume S''°'{k,uj)=M''°'{k,u)D{k,uj) 
with D{k, ei)=2/ {a+i-ei-i) [1] for calculations using this 
technique. We neglect instabilities in the ground state for 
simplicity and show intensities of DSFs in units of 1/J. 

This approach is applicable to the full range of mo- 
menta, energies, and magnetic fields for all directions of 
spin polarization except small energy scales O.IJ) and 
large \J'\ (> 0.7 J). The wide J' range is due to frustra- 
tion which prevents instability at ka:=7T, as will be dis- 
cussed later. The approach reproduces the 2D magnon 
above Hs and zero-field properties obtained in Ref. |10j . 

As a general property of the approach, spectral weights 
shift to higher (lower) energies for J'(fc)>0 (J'(fc)<0) 
[Eq. ID)] , and modes with 5-functional line shapes appear 
above (below) continua as poles of x""(fc, oj) in the RPA, 
which can be regarded as ABS (BS) of ID QPs 
A typical behavior is shown in the inset of Fig. [2fb). 
Although the modes get small intrinsic width due to very 
small weights around continua, the almost (5-functional 
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FIG. 2: S^\k,Lo), S+-{k,uj)/2, S" + (fc,a;)/2, and S{k,uj) 
(from above) at fcj,=0 for J'/J=0.34. (a) Results using the 
RPA at m=l/8, 1/4 and 3/8 (from the left), broadened in a 
Lorentzian form with full width at half maximum (FWHM) 
0.08 J. (b) Results using the method of Ref. [lOj] at m=l/4. 
Solid lines above (below) continua denote antibound (bound) 
states for Q'^k^Kir {TT<kx<2n). The inset shows the line 
shape at kx=1.25n. (c) Single-magnon modes in linear spin- 
wave theory at m=l/4. Here, the modes from two cone 
states, one with the ordering momentum +Q and the other 
with — Q, are shown, broadened in a Lorentzian form with 
FWHM=0.02J. 



peaks well represent QPs with integer S^. 

The results for DSFs using the RPA and the method 
of Ref. are shown in Figs. [DJa) and [2Kb), respec- 

tively. Strong intensities near upper (lower) edges of 
continua for J'(fc)>0 [0<kx<Tr] (J'(fc)<0 [7r<A:a;<27r]) in 
Fig. [2Ka) are signatures of ABS (BS), whose dispersion 
relations are asymmetric with respect to kx=T^ due to the 
sign change of J'{k) at kx=Tr. Each mode can be iden- 
tified as follows. In the top row of Fig. [D(a), the mode 
above (below) the continuum in S^^{k,uj) for J'(fc)>0 
( J'(fc)<0) can be identified as ABS (BS) of and ^^J* H], 
as indicated by solid lines in the top panel of Fig. [2Kb). 
[Dispersion relations as a function of ky for mod(fca;,7r)/0 
have a jump at J'(fc)=0.] For 5"' (fc,w), there are three 
kinds of continua as in the second row of Fig. ^a). The 
low-energy modes near kx=0 and 27r can be effectively 
regarded as QPs originating from the 1-0* mode in ID 
BMHIIll. The mode below the low -energy continuum 
near kx=Tr for J'(fc)<0 can be identified as BS of 2 ip*s 
The mode above (below) the high-energy continua 



near lu^2J can be effectively regarded as ABS (BS) of 
a and ip originating from 2-string solutions of the Bethe 
ansatz [i,[la]- For S '"(fc,a;), the mode below the con- 
tinuum for J'(fc)<0 can be identified as BS of 2 -^s 0] as 
in the third row of Fig. [2Ka). In the following, we dis- 
cuss distinctive features of these QPs. For comparison, 
behaviors of magnons in linear spin- wave theory 14, 23| 
are shown in Fig. ^c). 
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Multiparticle crossover.— We consider how QP pictures 
change from ID spinons interac ting through interchain 
exchange processes in zero field [10| to the 2D magnon 
above iJ^ In zero field, the mode of BS of spinons 
appears for J'(fc)<0 10], as in the top right panel of 
Fig. OJb). Yet considerable spectral weights remain 
spread in a broad continuum. Hence, the latter term 
in Eq. ([T]), which determines the ky dependence, has rel- 
atively small effects on the dispersion relation. Thus, 
signatures of ID spinons persist rather strongly in zero 
field In a magnetic field, the mode splits into three 
modes [BS of ip and tp* in S''''{k,u!), BS of a and V' 
in {k,uj), and BS of 2 V'S in S ^{k,uj)], as shown 
in the right panels in Fig. Eljb) and the left panels of 
Fig. [2l[a) near kx—l-5TT. As the magnetic field increases, 
these modes fade away as shown in Fig.[2Ia). Instead, the 
low-energy modes originating from the l^p* mode emerge 
near kx=0 and 27r in 5'"' (fc,a;), as in the second row of 
Fig. H^a). Using Eq. (P) and the l?/-* mode in ID, the 
dispersion relation is effectively expressed as 

cj(fc) = ei^.(fc,) H- J'(fc)M+o-(fc„ei^.)/2, (2) 

as shown by solid lines near kx=0 and 27r at low energies 
in the second panel of Fig. [Hb) . As the magnetic field 
increases, spectral weights concentrate in the modes, in- 
creasing Mj^(fc2;, eiv* ) in Eq. As a result, the modes 
become dominant and highly dispersing in both kx and 
ky directions like 2D magnons. Evidently, as seen above, 
the dominant modes in low and high fields are different in 
origin in contrast to magnons in linear spin-wave theory. 
Noting H=Lu{k=0) and using Eq.(l2]), we obtain 



(3) 



where i?iD is the magnetic field in ID [24 



This expres- 
sion is equivalent to that obtained in Ref. [25||. Compari- 
son with experimental results on CS2CUCI4 [15| is shown 
in Fig. [TJc). Hereafter, we use Eq. ([3]) to relate m to H. 

We confirm that the modes of Eq. ([2]) actually re- 
duce to the 2D magnon above Hg. By noting that 
the excitation energy in the Heisenberg chain above 
Hs is given as ei^. (^2,)= J(cos fcx — l)+i?iD, Eq. ([2]) re- 
duces to the dispersion relation of singie-spin-fiipped 
states 0: w(A;)= J(cos /cx-l)-h2 J'(cos ^cos 
with Eq. ([31). Small deviations from experimental results 
on CS2CUCI4 [Fig. [TJd)] will be due to other effects like 
Dzyaloshinsky-Moriya (DM) or interlayer interactions in 
the material 

IC ordering.— We consider instabilities at IC momenta 
expected from the behaviors of the BS of 1/' and ^* in 
S^^ik^uj). In the Heisenberg chain in a magnetic field, 
transverse low-energy correlations at kx=T^ are always 
more dominant than longitudinal ones at IC momenta 
[2^ . However, in the presence of frustrated interchain 
interactions with J' {kx^TT,ky)—0, spectral weights at IC 
momenta for J'(fc)<0 shift to lower energies by inter- 
chain interactions while those at kx='!T remain unaffected 




FIG. 3: Comparisons of line shapes of S{k,uj) with experi- 
mental results on CS2CUCI4 at ky—Q for (a) tj=0.35 meV and 
(b) kx — 1.5n. (Left panels) Solid lines are present RPA results 
of S{k, to) broadened in accordance with the experimental en- 
ergy resolution [^. Arrows in (b) indicate peak positions of 
S""{k,u}), a=+, z, and — (from the left). Symbols denote 
experimental results for H\\c in Ref. fsf. Here, the intensi- 
ties are rescaled after subtracting backgrounds [^, and the 
energies are normalized by J=0.374 meV [^. (Right pan- 
els) S{k,ui) using the RPA at ky=0 corresponding to the left 
panels, broadened in a Lorentzian form with FWHM=0.08 J. 
Green dashed lines indicate the scan paths for the left panels. 



at least in the present approximation [Eq. ([T])]. Thus, 
instabilities are expected at the momenta where the low- 
energy spectral weight in S^^{k,uj) becomes the most 
dominant i.e., k — (±7r(l -I- 2m) -I- 47rZi, 477/2) and 
(±7r(l — 2m) -I- 477/3, 27r + 477/4), where k {i—1^4) are in- 
tegers for J'>0, which correspond to the gapless points 
of the QP for BS oiip and ip*. The momenta for m=l/16 
are shown by solid circles in Fig. [ijb). Actually, IC or- 
dering whose momentum shifts as kx—Tr{l+2m) at ky—0 
has been observed in CS2CUCI4 for i?~l— 2 T||c [21] as in 
Fig- (He). Similar behaviors have also been observed in 
Cs2CuBr4 just below the 1/3 plateau [2^. These behav- 



iors have remained unexplained by spin- wave theory 14 ] . 
In linear spin-wave theory based on classically stable cone 
states, ordering occurs at the momentum of the gapless 
point in transverse DSFs [Fig.[2fc)], which is less depen- 
dent on H or m [3] ■ Note that by higher-order effects of 
J' for if~0 ^2^ and by large transverse correlations near 
Hs M , IS iSSi] , other states can be more stabilized than 



the IC orders. In real materials, ground-state properties 
are also sensitive to DM interactions l^, 14, 15, 2^, 2^. 

The above argument for the IC ordering will be ap- 
plicable to more general spatially anisotropic frustrated 
antiferromagnets with J' {kx=TT,ky)=Q. (For frustrated 2- 
leg ladders, we pick up 2 ky points.) Thus, the QP picture 
would be related to the IC ordering in frustrated mag- 
nets predicted using numerical simulations [13, 31, [34] 
and field-theoretical analyses 25, 3^ and discussed using 
an RPA for w=0 [13] ■ Note that the mechanism of the IC 



ordering is distinct from that of Ising-like chains 2^, 3S 
Excitation specira.— Inelastic neutron scattering exper- 
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iments can probe a quantity proportional to S{k, uj). We 
expect behaviors as shown in the fourth row of Fig. [UJa). 
Comparisons with experimental results on CS2CUCI4 Q 
for a;=0.35 meV are shown in Fig. [3)Ja). The line shapes 
can be interpreted in terms of the present QPs. For ex- 
ample, the shoulder near kx=2.6TT in the third row of 
Fig. [Slja) is mainly due to BS of 2 tps in S ^{k,uj) and 
BS of 2 iIj*s in {k,uj)- The peak near /c^^^Stt comes 
from ipip* excitations in S^^{k, oj). The dip near kx=3.4:TT 
is due to the shift of spectral weights near lower edges 
of ipijj* and 2ip continua to higher energies for J'(fe)>0 
[fcx>37r]. The asymmetry with respect to kx=3T: is un- 
derstood as a result of the sign change of J'{k) at kx—in. 

As in the right panels of Fig. ^h), the mode of BS of 
spinous in zero field splits into three modes in a magnetic 
field. The left panels show comparisons with experimen- 
tal results Q at fc = (I.Stt, 0). The low-energy peak is 
understood as superposition of BS of 2 ^ps in S ^{k,uj) 
and BS of V' and ip* in S^^{k, ui). Interestingly, the high- 
energy peak turned out to be a signature of BS of a and 
ip in (fc, uj) originating from 2-string solutions of the 
Bethe ansatz [oj, [16], which is not accounted for by the 
dominant modes in linear spin- wave theory [Fig. [UJc)] . 

In summary, through weak-interchain-coupling analy- 
ses using exact solutions of a chain, we introduced QPs 
for anisotropic-2D frustrated Heisenberg antiferromag- 
nets in a magnetic field, which can be regarded as BS 
or ABS of ID QPs. The QPs exhibit distinctive features 
different from those of ID QPs and those of magnons 
in linear spin-wave theory. Differences with ID QPs are 
obvious in view of their (5-functional line shapes with inte- 
ger S^, instability to IC ordering, and 2D features. Their 
statistics may also be different from those of ID QPs. Dif- 
ferences with magnons are the following. (1) The QPs 
show multiparticle crossover in a magnetic field. The 
dominant QPs in low fields fade away in higher fields. In- 
stead, other QPs which reduce to the 2D magnon above 
Hs emerge in a magnetic field. In linear spin-wave the- 
ory, magnons in low and high fields are the same in origin. 
(2) The QP in S^^{k, ui) causes instability to IC ordering 
whose momentum shifts as fc^=7r±27rm for J'(fc)<0 as 
in Fig. mb). In spin- wave theory, it is known that order- 
ing momenta do not shift so significantly as a function 
of iJ or m (3) The high-energy QPs in {k,Lu) 

come from 2-string solutions of the Bethe ansatz, whose 
behaviors are not explained by linear spin-wave theory. 

These novel features revealed in a magnetic field come 
from the stabilization mechanism of the QPs, which is not 
based on spontaneous symmetry breaking: in contrast to 
magnons created from classical spin configurations, the 
QPs are induced by interchain exchange processes from 
liquids of ID QPs. Thus, in analogy to collective modes 
in Fermi liquids, the QPs may be regarded as those in a 
kind of amsotropic-2D spin liquid or ip^ ip* ^ and a liquid, 
distinguished from the magnons. The arguments in this 
Letter can be generalized to weakly coupled antiferro- 



magnetic chains with J'{kx=TT,ky)=Q including ladders. 

Their relevance to real materials was confirmed 
through comparisons with experimental results. Various 
aspects observed in CS2CUCI4 in a magnetic field, such 
as the 2D dispersion relation, IC ordering, line shapes of 
S{k,uj), and the magnetization curve, were explained in 
a unified manner in terms of the QPs. Note that compar- 
isons in this Letter have no adjustable parameters except 
a single normalization factor in Fig. [31 It would be inter- 
esting in future studies to examine whether signatures of 
the QPs can persist in less frustrated or nearly spatially 
isotropic 2D antifcrromagncts in a magnetic field. 
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